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Abstract 



A minimal renormalizable SUSY SO(IO) model with B-L symmetry broken 
by 126 Higgs field has recently been shown to predict all neutrino mixings 
and the ratio Amg/Am^ in agreement with observations. Unlike models 
where B-L is broken by 16 Higgs, this model guarantees automatic R-parity 



(N 

, conservation and hence a stable dark matter as well as the absence of dim=4 

' baryon violating operator without any additional symmetry assumptions. In 

this paper, we discuss the predictions of the model for proton decay induced 
at the GUT scale. We scan over the parameter space of the model allowed 
. by neutrino data and find upper bounds on the partial lifetime for the modes 

^ ; r(n ^ TfOp) = 2t{p n+u) < (5.7 - 13) x 10^2 yrs and r(n K'^v) < 

P^; 2.97 X 10^3 yrs for the average squark mass of a TeV and wino mass of 200 

Q^' GeV, when the parameters satisfy the present lower limits on r(p — > K~^iy) 

^ ■ mode. These results can be used to test the model. 



I. INTRODUCTION 

Grand unified theories [1] provide an attractive framework for understanding the origin 
of the diverse strengths of the various forces observed in Nature. The basic idea is to have a 
single force associated with a grand unified local symmetry at a high scale, which below the 
scale of the symmetry breaking evolves into three different strengths corresponding to the 
observed weak, electromagnetic and strong interactions. The challenge is to have a theory 
where the three couplings evolved down to the Z-boson mass scale match their experimentally 
observed values. A concrete realization of this hypothesis is provided within the framework 
of supersymmetric models where the unification scale is Mjj ^ 2 x 10^^ GeV, provided one 
assumes that there are no intermediate scales below Mu and the scale of supersymmetry 
breaking is at the weak scale [2]. One may then assume that above the scale Mu, the grand 
unified symmetry takes over. While this observation suggests that the idea grand unification 
may be correct, one cannot take it as an evidence for it due to the extra assumption of grand 
desert which is crucial to obtaining unification of couplings. 

Typical grand unified theories not only unify forces but also quarks and leptons and 
thereby provide a hope to resolve some of the puzzles of the standard model. Especially 
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after the discovery of neutrino masses and mixings, there are a large number of observables 
in the matter sector that arc known fairly precisely and one may expect that a grand unified 
theory can provide an explanation of these observables. 

The obvious question then is: how can one test the idea of grand unification and more 
precisely, the existence of a local grand unifying symmetry of all matter and forces ? It is 
well known that one of the consequences of most grand unified theories is that proton is 
unstable and therefore one may argue that the true test of the grand unification hypothesis 
will come from the observation of proton decay. This has led to many experimental as well as 
theoretical investigations of proton decay in the context of simple grand unified models such 
as SU(5) and SO(IO). Although no evidence for proton decay has been found to date, the 
stringent experimental upper limits on the partial lifetimes to various modes have provided 
new and useful constraints on the nature of grand unified theories. 

The simplest grand unified theory is the minimal SUSY SU(5) model. In this model, 
the dominant decay of proton occurs via dimension five operators involving color triplet 
Higgsino exchange leading to the dominant decay mode [3] p — > K^v. The predictions of 
the minimal renormalizable SU(5) model for this mode has been discussed in many papers 
[5,6]^. The present experimental lower limit on this mode [9] is 1.9 x 10^^ yrs, which is 
an order of magnitude larger than prediction of the minimal renormalizable SUSY SU(5) 
model. Therefore this model is ruled out. It has been shown [6] that if one includes non- 
renormalizable terms in the superpotential [7], one can get somewhat higher lifetimes for 
this decay mode and the SU(5) model can still be consistent with experiments^. 

With the discovery of neutrino masses and mixings [10], SO(IO) is a much more inter- 
esting candidate for grand unification. It incorporates the right handed neutrino needed 
for implementing the seesaw mechanism [11] for understanding small neutrino masses. In 
addition, unlike SU(5), it incorporates all the fermions of the standard model as well as the 
right handed neutrino into one spinor representation. Furthermore, it is interesting that 
using the value of the neutrino mass difference required to understand atmospheric neutrino 
data and the seesaw formula, one can conclude that at least one of the right handed neutrino 
masses must be close to 10^^ GeV. The proximity of this scale to Mu suggests that they 
may indeed be one and the same. This would imply that seesaw mechanism may be the first 
signal of the idea of grand unification. It would then be urgent to look for other signals for 
SO(IO) grand unification such as proton decay [12]. 

It has recently been pointed out that there is a class of minimal SO(IO) models, where all 
neutrino masses and mixings can be predicted without assuming any additional symmetries 
[13-15]. The basic puzzle of neutrino physics i.e. the two large mixings (^0 < 9a) are 
explained in this model in a very interesting manner. All the Yukawa couplings that are 
responsible for proton decay are completely fixed in this model. Our goal in this paper is to 
study the proton decay in this model. 

The key ingredient of our model is the way B-L local symmetry in SO (10) is broken. 



By renormalizable, we mean a theory where only renormalizable terms are included in the 
superpotential. 

^For proton decay in string theories with SU(5) GUT, see [8]. 
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How B-L is broken decides whether the effective MSSM theory that emerges below the GUT 
scale conserves R-parity symmetry (defined hy (— or not. As is well-fcnown, 
R-parity symmetry is needed to guarantee the existence of a stable dark matter. The B-L 
symmetry could either be broken by the 16 Higgs or by the 126 Higgs multiplet. We are 
driven to using the 126- Higgs since, unhke breaking by 16, 126 breaks B-L by two units 
and leaves R-parity unbroken as can be concluded from the formula for R-parity given above 
and hence a stable dark matter [18]. 

The main observation in the papers [13-15] is that if B-L symmetry is broken by a 126 
dimensional representation of SO (10), the coupling of 126 to fermions unifies the flavor 
structure of the quarks and charged leptons with those of the neutrinos. This makes the 
model remarkably predictive in the neutrino sector and the predictions obtained in [15] are 
now in full accord with all known data for neutrino mixings. Furthermore, it is interesting 
that the experimentally much sought after parameter 1/^3 is predicted in the model to lie 
between 0.15-0.18, which can be probed in very near future in the MINOS experiment at 
Fermilab as well as several future experiments being proposed providing a test of the model. 

In order to study proton decay in this model, we consider all dimension five operators 
[16]. There are LLLL as well as RRRR type operators in this theory. We find that LLLL 
type operators dominate proton decay. This part of the discussion is similar to that of the 
SU(5) model. However, unlike the minimal SU(5), in the minimal SO(IO) model there are 
several contributions which for some domain of parameters can partially cancel each other. 
The cancellation is however not complete so that the net effect is to suppress the decay rate. 
We are then able to find upper limits on the proton lifetime. The lifetimes to various p-decay 
modes with charged lepton final states can be predicted as a function of the color triplet 
Higgsino mass and SUSY breaking parameters such as the wino mass and squark masses. 
For a specific choice of these parameters in the supersymmetry breaking sector (i.e. average 
Mq = 1 TeV and Mgiuino = 200 GeV), we find upper bounds for the partial lifetimes for 
the modes r(n ^ 7r°z/) = 2r(p ^ tt+z/) < {5.7 - 13) x 10^^ yrs and T{n ^ K^u) < 2.97x10^3 
yrs We also give the partial lifetimes for other charged lepton modes for these cases. 

This paper is organized as follows: in sec. 2, we review the basic outline of the model. 
In sec 3, we discuss the effective operators for the various proton decay modes, their origin 
and dependence on the parameters of the theory. In the first part of sec. 4, we present our 
predictions for the upper limits as well as the allowed ranges for the partial lifetimes to 
various modes. In sec. 5, we present our conclusions. 

II. BRIEF OVERVIEW OF THE MINIMAL SO(IO) MODEL WITH 126 

The SO(IO) model that we will work with in this paper has the following features: 
It contains three spinor 16-dim. superfields that contain the matter fields (denoted by 
ipa)', two Higgs fields, one in the 126-dim representation (denoted by A) that breaks the 
SU{2)ji X U{1)b-l symmetry down to U{1)y and another in the 10-dim representation (H) 
that breaks the SU{2)l x U{1)y down to U{l)em- These are the only two Higgs multiplets 
that couple to fermions and after symmetry breaking give rise to all the fermion masses 
including the neutrinos. The original SO(IO) symmetry can be broken down to the left-right 
group SU(2)l X SU(2)ji x f/(l) b-l by 54©210 Higgs fields denoted by S and S respectively. 
We wish to point out that the minimal realistic SO (10) model of the type under discussion 
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can be constructed without including a 54 Higgs field [17]. However we include it here since 
it provides the most general description of proton decay. 

To see what this model implies for fermion masses, let us explain how the MSSM doublets 
emerge and the consequent fermion mass sumrules they lead to. As noted, the 10 and 126 
contain two (2,2,1) and (2,2,15) submultiplets (under SU{2)l x SU{2)r x SU{^)c subgroup 
of SO(IO)). We denote the two pairs by 4)u,d and ^u.d^- At the GUT scale, by some doublet- 
triplet splitting mechanism these two pairs reduce to the MSSM Higgs pair (iif^, i^^), which 
can be expressed in terms of the and A as follows: 

if„ = cos q;„0„ + sin q;„A„ (1) 
Hd = cos ctd^d + sin aa^d 

The details of the doublet-triplet splitting mechanism that leads to the above equation are 
not relevant for what follows and we do not discuss it here. As in the case of MSSM, we will 
assume that the Higgs doublets H^^d have the vevs < >— vsinP and < >= v cos p. 

In orders to discuss fermion masses in this model, we start with the SO(IO) invariant 
superpotential giving the Yukawa couplings of the 16 dimensional matter spinor ipi (where 
i,j denote generations) with the Higgs fields Hiq = 10 and A = 126. 

Wy = hijipiipjHio + fijipiipjA (2) 

SO (10) invariance implies that h and / are symmetric matrices. We ignore the effects coming 
from the higher dimensional operators, as we mentioned earlier. Also we set all CP phases 
in the superpotential as well as vevs to zero, so that the observed kaon and B-CP violation 
is a consequence of all CP phases residing in the squark sector. 

Below the B-L breaking (seesaw) scale, we can write the superpotential terms for the 
charged fermion Yukawa couplings as: 

Wo = KQH^u'' + hdQHdd" + KLHde'' + txH^Hd (3) 

where 

h cos au + f sin (4) 
h cos ad + f sin ad 
h cos ad — 3/ sin ad 

In general ^ aa and this difference is responsible for nonzero CKM mixing angles. In 
terms of the GUT scale Yukawa couplings, one can write the fermion mass matrices (defined 
as jCjji — 'ipLMipji) at the seesaw scale as: 



^It must be pointed out that if the SO(IO) symmetry is broken by a 210 multiplet, then a 
new Higgs doublet pair from the (2,2,20)©(2,2,10) multiplets also mixes with the afore mentioned 
doublets. But this simply redefines the mixing angles a„,d and does not affect any of our results 



K = 
hd = 

he = 
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Ma = /iri + /r2 
Me = hr^ - 3r2/ 



(5) 



where 



h — h cos a 
/ 



sin Pvyjk 
f sin a„ sin Pvu,k 

COSftrf 



(6) 



^2 



COS au 
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sin a,. 



cot /3 
cot /3 



The mass sumrules in Eq. (6) were crucial to the predictivity of the model. The neutrino 
masses are then predicted using the type II seesaw formula [19]: 



m 



D 



(7) 



\Vb-L fVB-L 

In models which have asymptotic parity symmetry such as left-right or SO(IO) models, it 
is the type II seesaw that is more generic than the conventional seesaw formula. For some 
parameter range, the first term may dominate. What is important for our considerations 
is that the matrix / that determines the flavor structure of the neutrinos is related to the 
quark and lepton masses^. 

Using Eq.6 and 7, in reference [15] the neutrino masses and mixings have been calculated. 
We do not display those predictions here. But we note that by scanning over the allowed 
values for the extrapolated quark and lepton masses as well as quark mixings, we find a 
range of predictions for the neutrino sector. We find that a large range of the predictions 
are disfavored by the latest solar data [21]. However, there is also a significant allowed 
region. For this region, we extract all the Yukawa parameters hij and fij corresponding to 
this range and use them in our calculation of proton decay rate below. A typical set of 
values for /i's and fs in this range are: 

/3.26 X lO-'^ 1.50 X 10-^ 5.51 x 10-^\ 
-2.40 X 10-^ 
-0.0178 



h = 



1.50 x 10-^ 
V5.51 X 10-3 



5.51 X 10- 
-0.0178 
0.473 



(8) 



and 



/ 



/-7.04 X 10-^ 
-2.05 X 10"^ 
V-7.53 X 10-^ 



-2.05 X 10-5 
-1.85 X 10-3 
2.43 X 10-3 



-7.53 X 10-n 
2.43 X 10-3 
-1.64 X 10-3 y 



(9) 



^In extensions of the standard model with triplets, one has only the first term in the neutrino 
mass formula [20]. In these models however, the flavor structure of the neutrinos is completely 
independent of the quark sector. 
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FIG. 1. Superfield Feynman graphs that give rise to d = 5 effective proton decay operators. 



III. EFFECTIVE OPERATORS FOR PROTON DECAY 



In our model, there are four supersymmetric graphs that contribute to AB = 1 operator. 
They are given in Fig. 1 and involve the exchange of 10, 126 [22] Higgs multiplets and two 
mixed 10 — 126 diagrams. They will lead to both LLLL as well as RRRR type contributions 
given by the following effective superpotential: 



WaB=1 — Mj.^[Cijkieal3'yQ'iQjQkLl + 



r-) c,a ic,l3 c,7 d 



(10) 



where Mt is the effective mass of color triplet field. 

Note that one could in principle, diagonalize the mass matrix involving the color triplet 
superfields and write the Feynman diagrams in that basis. It is not hard to convince one self 
that the final result in this case will also have four parameters- an effective mass and three 
products of mixing angles. So by considering the above parametrization, we have not lost 
any information. This supersymmetric operator leads to effective dimension five operators 
that involve two quark (or quark-lepton) fields and two superpartner fields. In order for 
these operators to lead to a Four Fermi operator for proton decay, they must be "dressed" 
via the exchange of gluinos, winos, binos etc. Before we discuss this, let us first note that 
these operators must be antisymmetrized in flavor indices and then we get for the LLLL 
term 



^a(3-fMj.^[{Cijki — Ckjii)ufdju1ei — {Cijki — Cikji)ufdjd1h'i 



(11) 



There is a similar operator for the RRRR terms. As has been argued by various authors 
[4,23], for small to moderate tan/5 region of the supersymmetry parameter space, these 
contributions are smaller than the LLLL contributions. We also find this to be the case in 
our model. We will show this later; for the time being therefore, we will focus on the LLLL 
operator. 
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FIG. 2. Generic Feynman graph for dressing of d = 5 effective proton decay operators via 
gluino, Wino, Bino and Higgsinos. 

The effective four fermion operator responsible for proton decay can arise the gluino, 
bino and wino dressing of the above operators. The coefficient Cijki associated with the 
LLLL terms is expressible in terms of the products of the Yukawa couplings h and / which 
have already been determined by the neutrino and other fermion masses: 

Cijki = hijhki + xfijfki + yhijfki + zfijhki (12) 

where x, y, z are the ratios of the color triplet masses and mixings. As already noted, we do 
not need to know the detailed form for these parameters (x, z) in terms of these masses 
and mixings. In the end we will vary these parameters to get the maximum value for the 
partial lifetimes for the various decay modes. 

We now discuss the dressing of the various terms. The typical diagrams are shown in 
Fig.2. 

A. Gluino dressing 

It has been pointed out in several papers [24] in the limit of all squark masses being 
same as in mSUGRA type models, these contributions to the effective four-Fermi operator 
for proton decay vanishes. It results from the use of Fierz identity for two component spinors 
which states 

(0102) (0304) + (0103) (0204) + (0104) (0203) = (13) 

0j are the chiral two component spinors representing quarks and leptons and (A-B) = 
= e'^^A^Bp where a and /? are the spinor indices (a,/? = 1,2). Since satisfying the 
flavor changing neutral current (FCNC) constraints allow only very small deviations from 
universality of squark masses, the gluino diagrams should be small (proportional to ^LL.ii in 
standard notation [25]) in realistic models. We will therefore ignore these contributions. 
The same results hold also for the RRRR operators. 

B. Neutral Wino and Bino contribution 

To analyze the contribution from and S, we choose the the operator fig = UfDjU^Ei 
as an example. Note that we can use W° and B in the loop instead of the the superpartners 
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of Z boson and photon is because they are both mass eigenstates due the the assumption of 
the universal mass. 

1. B dressing 

There are 6 different dressings of the operator VL^ through B. We can spht them into 
two groups. One group involves the lepton and the other one does not. Within each group, 
the product the hyperchargcs from the two vertices are same. Each of these groups then 
gives zero due to the Fierz identity as in the case of gluino dressing. This show that the B 
dressing is zero by the same Fierzing argument as the gluino case in the limit of universal 
squark masses. 

2. W° dressing 

For the W° case, the vertex involving the lepton is same as that of quarks but different 
by a negative sign between up type and down type particles. The dressing of uu and dd are 
then different from that of ud by a negative sign. Because the W° are lepton/quark blind 
and the dressing does not change anything except from boson to fermion, the two groups 
we used in the B analysis are the same. So after dressing, we have 

- 2(-«df)KeO - «e,)K4) + (14) 

By the Fierz identity, the sum of the first two terms is equal to the third and so we have 

n,^A{uyl){d^^ei) (15) 

Due to the antisymmetry of this expression in the color indices, it is antisymmetric in 
the interchange of i and j. This implies that i must be different from k and so the two up 
quarks belong to different family. This antisymmetry remains true even after we pass to the 
mass eigenstate basis, as is easily checked. The result is simply due to {u^J.) — —{u'^uj). 
The conclusion is that there is no + or tt" + ef decay mode from the W° dressing. For 
the same analysis, the operator U^DjDlui gives A{dPjC[l){u'^'Ui) and so it only contributes to 
+ vi decay mode. 

C. Wino contribution 

In view of the discussion just given the dominant contribution to proton decay arises 
from charged wino exchange converting the two sfermions to fermions. These diagrams have 
been evaluated in earlier works [5,6]; we will assume that all scalar superpartners have the 
same mass. This leads to the following effective Hamiltonian: 

>Cab=i 2Ieaf3^{Ckjii- Cijki)[ul''^Cd'^d]''^Ciyi + uf^Cdlu"''^Cei\ (16) 

where / is given by / = Using this expression and adding a similar contribution 

7 

from exchange, we can now write down the C coefficients for the different proton decay 
operators. Table I lists the total contributions to the different operators in the leading order: 
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Table I 



Operator 


C-coefficient 




2Isme,(C2ni-Ci,2i) 




2/(Cii2; — C1211) 




2/sin6'c(C22i/ - C212Z) 


udue^ 


2/sin^e(C2iiz-Cii20 


usde£ 


2/(Cii2Z — Ci2iz) 



Table Caption: The coefficients for various AB = 1 operators from tlie GUT theory. 
The C's are products of the Yukawa couplings in the superpotential as in Eq. (12). 

D. Estimates of the RRRR operators 

In this subsection, we give an estimate of the RRRR operators and confirm that they 
are indeed negligible compared to to the LLLL operator contributions for moderate tan /3 
region that we are interested in. First we note that the gluino dressing graphs are zero in 
the limit of all squark and slepton masses being equal, by the same argument as for the 
LLLL operators. Secondly, since all superfields in this operator are SU (2) ^ singlets, there 
are no wino contribution to leading order. The only contributions are therefore from the 
bino exchange and the Higgsino exchange. 

Bino exchange generates a four Fermi operator of the form eaP'yuf''^C(ffJu1^''^Cei. (where 
c in the superscript stands for charge conjugate). This operator in the flavor basis must be 
antisymmetric in the exchange of the two flavor indices i and j. Once they are antisymmetric 
in the flavor basis, they have to involve charm quark in the mass basis since uu terms will 
then be zero. Thus the to leading order the bino contribution also vanishes. 

The Higgsino exchange leads to an effective operator of the form: 

Ieaj3'y{Dkji'i')Xi'i,i'i[ul'^''^Cdf{d]'*Cul) + uf^Cdlui^'^Cei] (17) 

where Xi'i^in ~ ^^^2^ ^ p Mu^i'iMg^in . Since 1/ sin/? cos/? ~ tan/? for large values of tan/?, 
this contribution grows with tan /?. It is clear from inspection that the largest value for this 
amplitude comes from t intermediate states and we estimate the largest contribution to be 
of order Cis23^^^^^j^ — 10~^° as compared to the LLLL contribution which are of order 

Cu23^ ~ 10~^. Therefore, we can ignore the RRRR contribution in our discussion. 

IV. PREDICTIONS FOR PROTON DECAY 

Let us first note that the operators with s quark lead to p-decay final states with K meson 
whereas the ones without s lead to tt final states. Also generally speaking the amplitude 
for nonstrange final states are down by a factor of Cabibbo angle (~ 0.22) compared to 
the strange final states as in the case of SU(5) model. However, as we will see, we need to 
do a fine tuning among the parameters x, y, z to make the p — > + 1? compatible with 
experiments. The same fine tuning however does not simultaneously lower the amplitudes 
with nonstrange final states. As a result for some domain of the allowed parameter space. 
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one can have the p — > tt"*" + P mode as the dominant mode. This is very different from the 
minimal SU(5) case. 

In order to proceed to the calculation of proton lifetime, we must extrapolate the above 
operators defined at the GUT scale first to the Mg and then to the one GeV scale. These 
extrapolation factor have been calculated in the hterature for MSSM and we take these 
values. The required factors arc: AiAg [4] and are given numerically to be Al — 0.4 (SUSY 
to one GcV scale) and Ag = 0.9 - 1.0 (GUT to Ms scale). 

The next step in the calculation is to go from three quarks to proton. The parameter 
is denoted in the hterature by (3 and has units of (GeV)^. This has been calculated using 
lattice as well as other methods and the number appears to be: /3 ~ 0.007 — 0.028 [26]. We 
find that for our choice of the average superpartner masses, for f3 > 0.01, there is no range 
for the parameters {x, y, z) where all decay modes have lifetimes above the present lower 
limits. Of course as the superpartner masses increase, larger j3 values become acceptable. 
For instance, we note that a change Smym^ by 10% allows a 20% higher value in /3. We 
confine ourselves to the domain 0.007 < (P/GeV^) < 0.01 and find that for all choices of the 
free parameters allowed by the present lower limits, lifetimes for the decay modes p — > vr"'"^ 
and n — > tt'^P have upper limits, which can therefore be used to test the model (see below). 

Finally, in a detailed evaluation of proton decay rate to different final states, we take 
into account the chiral symmetry breaking effects following a chiral Lagrangian model (the 
first two papers of Ref. [26]), where the chiral symmetry breaking effects are parameterized 
by two parameters D and F. These are usually chosen to be the same as the analogous 
parameters in weak semileptonic decays [29]. 

For this case, we find the rate for proton decay to a particular decay mode P£ [P is the 
meson and ^ denotes the Icpton) to have the form: 

where /(F, D) is a factor that depends on the hadronic parameters F and D and we have used 
13 — 0.01 GeV^ in the last expression. We now discuss the evaluation of the parameter |Cp 
which determines the partial proton decay lifetimes for various modes. The relevant modes 
are p K^u, K^fi^, K^e~^, Tre"*", tt/j,^. The present lower limits (including n —>■ ttv^Ku 
modes) on these modes are: 



Table II 



mode 


hfetimes (xlO^^ yrs) 


p — >■ K^v 


19 


p K'>e+ 


5.4 


p K^n+ 


10 


P — >• TT+P 


0.2 [27], 0.16 [28] 


p TT^e'^ 


50 


p —f n''^fi^ 


37 


n TT^U 


4.4 




1.8 
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Table caption: Present experimental lower limits on the relevant proton decay modes 
from Super-Kamiokande and Kamiokande experiments. 

To proceed with this discussion, first note that C's are products of the known Yukawa 
couphng parameters h and / and the four GUT scale parameters as already discussed in 
Eq.(12). The GUT scale values of h and / are obtained from neutrino fits described in Ref 
[15] and are given in the previous section. 

As far as the GUT scale parameters go, we will keep the overall mass parameter to be the 
GUT scale i.e. 2 x 10^^ GeV. We have diagonalized the mass matrix of the color triplet GUT 
scale Higgs fields in 10, 126 etc and we find that they also lead to the same parametrization 
as we have given here. The meaning of the overall mass scale is then that it represents a 
product of one of the mass eigenstates with the determinant. We have checked that for the 
allowed range of parameters, the value of the determinant, given by \x — yz\ is around 0.25 
or so, so that none of the mass eigenstates is too much higher than the GUT scale. As a 
result the threshold effects on the gauge coupling unification is minimal. 

We then adopt the strategy that we vary the parameters x,y,z in such a way that the 
nucleon decay rate to the p K^v mode (summed over all the final neutrino final states) is 
consistent with the present experimental lower limit. Since there are three final states which 
add incoherently, this narrows the space of the x, y,z io a. small domain. In this domain we 
pick a point (call it (xq, ^/o, Zq)), where all other modes also satisfy their present experimental 
constraints as in Table II. We then vary the (x, y, z) parameters around (xq, yoi Zq) until the 
lifetime for a mode goes below its present experimental lower limit. We find that dependence 
on the parameter z is much stronger than the others. In Fig. 3 and 4, we give the allowed 
domain of the parameters (a;, y) consistent with the various experimental lower limits on the 
partial lifetimes for an optimum value of z. The boundary of the domain is determined by 
the lower limit on the the p — > K^V. Inside this domain the t{p ^ K^i>) is higher than 
its present lower limit. The maximum value of the p — > 7r"'"P and n tt^u occurs at the 
boundary. We find that r(n n^i?) = 2r(p tt+z/) has an upper bound of (5.7—13) x 10^^ 
yrs depending on whether /5 = 0.01 — 0.007 GeV^. At a different point in the parameter 
space, r(n Kv) acquires its maximum value of 2.9 x 10^^ years. The predictions for the 
partial hfetimes of other modes are given in Table III for both these cases. These values are 
accessible to the next round of proton decay searches. 



Table III 



mode 


r/1032 ^ ^ 0.01 


r/1032 yrs: {3 = 0.007 


r/1032 yrs: (3 = 0.007 




r(n — > ttP) maximized 


T(n — > ttP) maximized 


T(n — > Ku) maximized 


p K^v 


19 


19 


19 


p K'^e+ 


1793 


2848 


188 


p K^fi^ 


184 


303 


28 


P — >• Tl'^U 


2.87 


6.5 


2.59 


n — > 7r°z/ 


5.7 


13 


5.18 


p — > 7r°e+ 


2452 


3857 


243 


p — > 7r°/i+ 


263 


430 


37 




1.9 


3.1 


29.7 



Table caption: Predictions for various nucleon decay modes for the case when the lifetime 
for the mode n ^ 7r° + P attains its maximum value. The units for (3 parameter (i.e. 
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GcV"') has been omitted in the table. In column 4, we give the lifetimes for the case when 
r(n Kv) is maximized. 

We check the above results adopting an alternative strategy where we express the three 
parameters [x, y, z) in terms of three partial hfe times and plot the other lifetimes as a 
function of these partial life times. It turns out that if we pick a certain value for the partial 
life time of the p — > K^x mode and use it as an input, the other two input values get very 
restricted. This allows us to use only the p K^fi^ mode as a variable and give the others 
as a prediction. In Fig. 5 and 6 we present the allowed values for various partial lifetimes as 
a function of the partial lifetime for the mode p — > K^^^ . There is a slight spread around 
the various lines. We first find that the lifetime for the mode K^V can be arbitrarily large 
as can be seen from Fig. 5. Also, from Fig. 5, we see that modes n — > 7r°P and n — > K^D 
have upper bounds which are same as the ones derived previously. 

V. CONCLUSION 

In summary, we have discussed the predictions for nucleon decay in the minimal SO (10) 
model which has recently been shown to lead to predictions for neutrino masses and mixings 
in agreement with observations using only the assumption of type II seesaw mechanism. 
The key feature of this model is that B-L symmetry is broken by a single 126 field that 
also contributes to fermion masses. For the range of the parameters that are allowed by 
the neutrino data, we vary the GUT scale parameters ( unrelated to the neutrino sector) so 
as to satisfy the stringent experimental bounds for the decay mode p — > + P. We then 
predict an upper limit for the lifetimes for the modes p ^ -k^ + v and n — > 7r°P as follows: 
T{n 7r°p) = 2r(p tt+z/) < 5.7 — 13 x 10^^ years and r(n Kv) < 2.9 x 10^^ yrs for 
the wino masses of 200 GeV and squark and slepton masses of a TeV. This should provide 
new motivations for a new search for proton decay, more specifically for these decay modes 
in question. 
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FIG. 3. Allowed Region for (x, y) coming from experimental lower limits on lifetimes for differ- 
ent decay modes for z = 0.329. The point (Aa;, Ay) = (0, 0) corresponds to (x, y) = (—0.036, 0.387). 
Note that the region is most constrained hy p ^ K + v' mode. 
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FIG. 4. Upper limit on the n ^ n + v partial lifetime while satisfying bounds on the lifetimes 
of all other modes. The point (Ax, Ay) = (0,0) corresponds to {x,y,z) = (—0.132,0.347,0.306). 
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-2 □ 2 H. ife 

Log(i:K^/r^) 

FIG. 5. This figure gives the values of the lifetimes for different proton decay modes as a 
function of the lifetime of thep — >■ K^fi'^ mode (represented here by logio^-^ where tq = 14.6 x 10^^ 
years) when r(p K^v) mode is at its maximum value. This figure displays the values for one 
range of (x, y, z) and the following figure does it for a complementary range. Also note that we 
have not included the gauge boson exchange diagrams, which provide a value for these lifetimes 
around 10^^ years or so. 
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